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Explicit chiral symmetry breaking is a natural feature of many QCD inspired models with multi- 
quark interactions. To carry out the 1/N expansion of these theories, one integrates over the quark 
fields. In the process of calculating the one-loop fermion determinant which is not chiral invariant 
C | due to quark masses, special care must be taken in order not to alter the original symmetry breaking 

, pattern of the Lagrangian. We show how to do this consistently. 
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The purpose of this paper is to resolve a problem that arises in models with chirally symmetric multi-fermion 
interactions when this symmetry is explicitly broken by the mass term of the fundamental fermion field. These can 
be, in particular, the corresponding extensions of the SU(N) Thirring model in two space-time dimensions [H, or 
the Nambu - Jona-Lasinio like models in four space-time dimensions 0, [I[ . To carry out the 1 /N expansion of these 
theories, one should integrate over the fermions, facing the task of calculating the one-loop fermion diagrams to obtain 
the effective action of the auxiliary bosonic fields. The problem is that the standard procedure leads to a result which 
does not possess the transformation properties of the original fermion Lagrangian with respect to the action of the 
continuous chiral group. 

The essence of the problem has been already discussed in Ref. |4| , where the formal solution has been found by 
adding a functional P in the bosonic fields and their derivatives to the real part of the fermion determinant. Here we 
suggest an alternative method to solve the problem. This new solution leads directly to the result and in addition 
1^ ' provides for an interpretation of the correcting functional P in terms of Feynman diagrams. 
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I. INTRODUCTION 



II. THE 1/N EXPANSION OF THE SU(N) THIRRING MODEL WITH EXPLICITLY BROKEN 

CHIRAL SYMMETRY 



OO 

o 

The SU (N) Thirring model describes a system of N Dirac fields tpr a -j carrying the flavor index a — 1,2,... JV, which 
we promptly suppress, in two space-time dimensions with the Lagrangian density 



C(yj, rl) = _^(<0 + m)i> + 9 — [(^) 2 + (^V') 2 ] , (1) 

where ip — $ — 7 M <9 M . The 2x2 antihermitian symmetric matricies a M = —ifij^, [i = 0, 1 generate the Clifford 
algebra: {a^,a u } = 2<5„„. The matrix (3 is imaginary and antisymmetric. We take here the choice j3 — — 7 = u 2 , 
7 1 = iai, 75 = 03, where ui are the Pauli matricies. One has also {7^,7"} = 2g^ v with the Minkowski metric tensor 
€T = diag(l, -1). 

The Lagrangian possesses a ?7(1) chiral symmetry ip/ a \ — > e ie75 i/>( a ) which is explicitly broken by the mass term. 
One can easily find that under small variations 9 <C 1, we have 

6 , %p) = —fhd^ipip) = —26rh'4>i r )5il). (2) 

Following Ref. [H, one can equivalently introduce auxiliary fields a, <fr, and write 

C(lp, tp, cr, (p) = -4>{i$ + m + a + 475^)^1 - - . (3) 



The Euler - Lagrange equations for the bosonic fields are constraints 

-?£ = fah + — = to -— 

da g 2 d(p 



— = 0, -— =^75^ + 4 = °! ( 4 ) 
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which relate the chiral transformations of fermions and bosons, and hence 



6a = -<? 2 <W) = 200, (5) 
6(f) = - g 2 5(ipi l5 ip) = -26a. (6) 



Then it follows from ([2j) and ([5|) that the infinitesimal symmetry transformation is now 



6C=^ 6a. (7) 
9 

Performing the Gaussian integral over the fermions in the functional integral corresponding to §5§ (see Eq. (fT5|) 
below), one finds that the real part of the effective action S e f{ has a different transformation property. Indeed, one 
obtains in euclidean space 

Re Sf = ^ Tr In D E D E -J d 2 x E ^-f- (8) 

with 

6(ReSf)=mNTr(-^-), (9) 
\D E U E ) 

where the Dirac operator is given by D E — — (ift + fh + a + ij5(f>). Here we use the following convention: a Lorentz 



-ix 



xf . The euclidean 7-matrices are antihermitian 7 — » 7^ = — i&2, 
7 1 — * -ff = ia\. They satisfy the anticommutation relation {7^,7^} = — 25 a b- Accordingly ft — j^d^ is a hermitian 
operator. Furthemore Tr = fd 2 x E tr, where the last trace is over Dirac indices. 

The unsuccessful outcome of Eq. © shows that the standard evaluation of the Gaussian functional integral over 
fermions, when chiral symmetry is explicitly broken, must be corrected, otherwise the result is not consistent with the 
constraint imposed by the symmetry breaking pattern of the theory at the tree level. In particular, this can be done 
by the insertion of new counterterms into the effective action ijHJ). The corresponding method has been developed 
in However, it was not clear then (a) if these counterterms were compatible with the dynamics of the original 
multi-fermion system, and (b) if the method could be applied to renormalizable models. In the following we shall find 
the answers to these questions describing a new way of dealing with the problem. 



III. TADPOLE MECHANISM 



It is convenient to divide the Lagrangian ([3]) into two parts C = C\ + £2 ■ The first part 

_ TYICT 

£1 = -iptift + rh + a + i^M 5- (10) 

9 

is invariant under the action of the chiral group. We expect that this property will be still fulfilled for the corresponding 
part of the effective Lagrangian obtained as a result of integration over the i/j fields. Our expectation is based on the 
observation that the dangerous symmetry breaking fermion mass term can be easily subtracted if one considers the 
tree-level tadpole a — > vacuum transition contained in Cx, as it is shown in Fig. [1] 




FIG. 1: The mass term of the fermion is cancelled due to a contribution generated by the a — » vacuum transition. 
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The second part 
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(11) 



possesses the necessary transformation property 8C2 = rh5a/g 2 of the Lagrangian Note, that the term fha/g 2 , 
which we subtract in (jlOp and add in is unambiguously determined by the chiral symmetry restriction 8C\ = 0. 
It follows from Li additionally that the auxiliary fields a and cf> do not propagate at tree level. They have the 5-like 
Green functions 



-v) = -ig 2 6(x - y). 



(12) 



To integrate the fermions out of the Lagrangian L\ and not break its transformation properties with respect to the 
chiral group, we modify the corresponding functional integral, writing it as 



VipVip < exp 




A <7 (x 1 - x 2 ] 



ST(x 2 ) 



exp ( —i / d 2 x[ip(i$ + fh + a + 175^)^ H — ~T] 



|r=o 



(13) 

Obviously, this definition coincides with the standard formula at fh = 0. For fh ^ the action of the first exponent 
upon the second is to include as basic element in the functional integral the chiral symmetric tree-level result shown 
in Fig. [TJ On the other hand, by first integrating out the fermions, this exponent generates an infinite number 
of all possible one-loop fermion diagrams attached to the external auxiliary fields, including those that contain the 
a — > vacuum transitions, rendering the result invariant under chiral transformations. Some of these diagrams are 
shown in Fig. O 
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FIG. 2: Diagrams which result when the first exponential operates on the functional integral (|13[l . Full lines represent fermion 
propagators, dashed lines with a cross a — > vacuum transitions, dashed lines stand for either a or <f) fields. 
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Indeed, taking into account eq. (fT2|) . one easily integrates over X2 in (fl3l) . Then, using the formula 
5 



st( Xi) j exp {-'J d x 7 r{x) ) = ) cxp rv d * ? r(x) J ' (14) 

we derive for (|13[) an expression, where after setting T = and taking the first exponent out of the functional integral 
we get 



exp y-fh J d 2 x - - j ) J VipVip exp y-i J d 2 x[ip(i$ + m + a + ij 5 4>)ip] J . (15) 

Since the first exponent is a translation operator, which shifts the argument a — » cr — m, we obtain finally for (|13[) 

V$V^exp(-ijd 2 x[ip(i@ + <T + ij 5 (j>)i>]J <x(detD E \~ =Q ) N . (16) 

Thus, we arrive at the following bare effective action 

**Sf = f Trlnrt E D E \^ -/d^ ^|^ , (17) 

where the real part of the fermion determinant is invariant with respect to the chiral transformations ([5])-©, as we 
wished to find, and SC c s — fh6a/g 2 , as it follows from the second term. 

One might think of criticizing the above calculations by claiming that the same result (that is Eq. (JT7J)) can be 
obtained by the usual method, i.e., by shifting a — > a — fh in ([3]) and subsequently integrating over the fermions. 
We argue, however, that the shift of the scalar field in the generating functional, being a replacement of the variable, 
changes automatically the transformation law of the pseudoscalar <fi to 8<j) = —28(<t — fh). As a consequence, the 
symmetry breaking pattern of the effective Lagrangian obtained in this way differs from |(7j) and, therefore, suffers 
from the problem indicated in Sec. II. In contrast, Eq. (TIT)) obtained through the complementary tadpole diagrams 
contained in the definition Q13[) leaves the auxiliary fields and their trasformation properties unchanged. 

IV. TADPOLES AT THE LEVEL OF THE GENERATING FUNCTIONAL 

The procedure considered above can be easily implemented in the framework of the functional integral approach. 
Indeed, writing the linearized Lagrangian density C(ip,ip,cr,<f>) as the sum of the symmetry breaking, Cs-B^^ip), and 
symmetric, £si parts, one has for the generating functional 

Z = j Thplhi)expiS(<il},i>) = [ V^V^ e lSsu(i [ VaV<l> exp iS s (ip,$, <?,<!>), (18) 



where 5sb(V':V') = — J d 2 x ipfhip . At this stage, to carry out the l/N expansion, one usually integrates over the 
fermions assuming that the order of integration in (|18p can be changed without serious consequences. This is not 
true, actually. The reason is that the double integral over bosonic fields is a chiral invariant, but the one over fermions 
is not. Thus, the pattern of symmetry breaking will be altered. 

We argue now that the consistent procedure here should be as follows 

Z = J VaVcj)e iSaii(a) J T>ipVi>expiS s (i},$,a,(t>), (19) 

i.e., the change of the order of integrations in Ijl8|) should be accompanied with the corresponding replacement 
of the symmetry breaking term: Csb^j w) ~^ — 7 '9' '■> where £sb(c) is the solution of the equation 

<5£sb(o") = <5^sb(V': V0- This equation can be solved, because chiral transformation properties of bosonic and fermionic 
fields are mutually correlated. As a result the functional integral over fermions in Eq. (|19[) is chirally invariant 
(excluding the anomaly, which is not important for the question studied here), and therefore the symmetry breaking 
pattern is strictly traced. 

At first sight, Eq. (| 19[) looks as if we did the shift cr — * a — fh in (fl"8|) to obtain it. This is not the case. As we have 
already learned, the shift does not change the symmetry breaking pattern of the fermionic part of the Lagrangian, 
and therefore cannot lead to the desired chirally symmetric functional integral over fermions. Eq. (|19|) should instead 
be considered either as a postulate related with the symmetry or as the result of the c-tadpole mechanism. The latter 
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is in our understanding the reason for the unconventional rule associated with the order of integrations in Z . This 
equation contains the main message of our work because it expresses the result of Sect. IIIII in the shortest and most 
general form. 

One can recognize in the replacement £sb(V , :V') ~ * ^Sb(c) the use of the constraining relation ((?]). On one hand, 
this explains why Eq. (fl~9]l is inwardly consistent with (fT8|) . On the other hand, if one accepts that it is this field 
equation that controls the step from (I18|) to (fT9|) , the whole procedure can be straightforwardly applied to any theory 
with broken discrete 75-symmetry. The massive GN model [5| is an example of such a theory. In this particular case 
our formula leads to the result obtained by Feinberg and Zee for the effective potential. 

Next, let us try to understand in simple terms what is the main difference between our calculations here and the 
method presented in Following those works, we would integrate out the fermions directly in Eq. (fT5)l . Of course, 
we would not pay much attention to the order of integrations there. However, we would calculate systematically (in 
the framework of the heat kernel expansion) the counterterms which are needed for consistency with the requirements 
of chiral symmetry and would add them to the effective Lagrangian. What would be the result? The answer is very 
simple. It would coincide with the effective Lagrangian which one finds integrating out the fermions in Eq. (|19| . 
To see this it is instructive to look again at Eq. (fl~3|) . where the first exponent, being taken out of the integral and 
applied to the result of integration, generates the infinite set of counterterms which are necessary to recover the 
correct symmetry breaking pattern of the outcome. It is exactly this what has been done in [4j, but with the use of 
a different and more complicated technique. It is clear that both results coincide. 



V. FINAL COMMENTS AND CONCLUSION 

The renormalization of the effective Lagrangian (jTTJ) can be done similarly to the massive GN model In order 
that chiral symmetry constraints be valid in the renormalizable theory, it is necessary that all divergences of the 
theory may be absorbed in two constants, the mass fh and the coupling constant g 2 . It is also important that the 
boson variables are chosen such that their chiral transformations do not depend on m, otherwise one would have 
fatal infinities in the symmetry transformation laws. These conditions are fulfilled for (|17p and we conclude that our 
calculations do not destroy the renormalizability of the multi-fermion Lagrangian in two space-time dimensions. 

The method presented here can be easily extended to the known 4-dimensional models with nonabelian chiral 
symmetry and multi-fermion interactions. Several examples are given by the Nambu - Jona-Lasinio type models of 
QCD with explicit chiral symmetry breaking. In these cases the original quark Lagrangian satisfies the hypothesis 
of partially conserved axial vector current (PCAC) at the quark level and the tadpole mechanism which we suggest 
protects the theory from loosing this property in the large N c limit. There are several examples in the literature where 
the subtle but essential difference between Eqs. (|18[) and (fH?)) was not recognized, leading to an unjustified violation 
of the PCAC relation. Note also that the lattice calculations for the massive GN model show that PCAC is satisfied 
in the continuum limit @ . 

To conclude, we have outlined the basic steps required to formulate the functional integration over multi-fermion 
interactions with explicit chiral symmetry breaking in leading 1/iV order. Without loss of generality we have considered 
the massive GN model, which is being actively studied nowadays • We have shown that the method preserves 

the symmetry breaking pattern of the original Lagrangian and does not alter its dynamics. The information over all 
counterterms which are needed to be added to fulfill these requirements is cast into a very compact form and readily 
implemented as a rule in the functional integration. 
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